PROCEEDINGS 

OF SCIENCE 



o 
o 

(N 

Or 
<D 

(N 



Dynamical simulations with HYP-link Wilson 
fermions 



Oh! 



> 

o 

m 

O 

o 



Stefan Schaefer* 

NIC, DESY, Platanenallee 6, D-15738 Zeuthen, Germany 
E-mail: 



Stefan . schaef er@desy . de 



Anna Hasenfratz, Roland Hoffmann 

Department of Physics, University of Colorado, Boulder, CO 80309 USA 



E-mail: anna@pizero . Colorado . edu hof fmann@pizero . Colorado . edu 



We present results from simulations of two dynamical flavors of improved Wilson fermions with 
nHYP smeared gauge links. We demonstrate that the simulation is stable at a pseudo-scalar mass 
of 360MeV, a 2. lfm box and a lattice spacing of 0. 13fm. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Dynamical simulations with HYP -link Wilson fermions 



Stefan Schaefer 



1. Introduction 

Ever since the introduction of the original APE smearing in 1987[jl|], many applications have 
shown the beneficial effects of constructing the lattice Dirac operator from smeared links. The 
hypothesis behind employing smeared links is that short range fluctuations of the gauge field, the 
so-called dislocations, are responsible for some of the poor scaling behavior and large cost of 
fermion simulations. For different lattice discretizations of the Dirac operator, the improvement due 
to smearing comes in different disguises. In Wilson fermions dislocations cause the exceptionally 
small eigenvalues of the Dirac operator, in staggered fermions the large taste breaking. Domain 
Wall fermions acquire their residual explicit chiral symmetry breaking through dislocations and 
they are responsible for the high numerical cost of constructing the overlap operator. 

An alternative approach to reduce these effects is the use of special gauge actions, e.g. the 
Iwasaki action or DBW2, which suppress the occurrence of dislocations. However, it has turned 
out that these gauge actions themselves induce quite large scaling violations. Smearing the links, 
however, reduces the effect of the dislocations on the fermions only. They remain part of the gauge 
dynamics and one can choose a gauge action which does not introduce poor scaling behavior from 
the gluonic sector. Increased auto-correlation times in molecular dynamics based algorithms have 
been observed when using improved gauge actions as well. Below, we will show that we do not 
observe signs of this in our simulations. 

Regardless of whether one accepts the explanation of these effects, in many quenched studies 
it was demonstrated that smearing helps to improve scaling in the situations listed above. It is 
therefore natural to use it in dynamical simulations too. 

The goal of this conference contribution is to convince the reader that simulations of improved 
Wilson fermions constructed from the recently suggested nHYP links [^] are stable even at a coarse 
lattice spacing. The additional cost of the smearing is small and more than compensated by the 
improved conditioning number of the fermion matrix. Exploratory studies of non-perturbative 
improvement of this action using Schrodinger functional techniques have also been presented at 
this conference [Q]. 



2. Smearing procedure 

Most dynamical algorithms are based on molecular dynamics and therefore one needs to dif- 
ferentiate the action with respect to the gauge fields. If one uses smeared links, the derivative of 
the smeared link with respect to the thin link is needed. This turns out to be a problem for the 
projection which is part of the definition of the APE smearing, 



V n ,n = Proj 



SU(3) 



±v^ 



(2.1) 



with n labeling the site and the result of the projection B = Proj SU ( 3 )A is defined as the matrix 
B E SU(3) which maximizes tr[A^B + B^A]. The stout smearing of Morningstar and Peardon [Q] is 
a fully differentiable alternative. However, we found it to be considerably less effective in reducing 
the effect of the dislocations. 
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It also turns out that typically one level of smearing is not enough and one needs to iterate the 
procedure. However, iterating it too many times leads to Dirac operators with a large footprint. 
In quenched studies, a particular smearing recipe has proved to be efficient but not over-doing it: 
HYP smearing^]. It consists of three levels of projected APE smearing, however, the smearing 
is restricted such that the smeared link only receives contributions from its hypercube. In Ref. [||] 
we introduced n-HYP smearing which differs from the original HYP smearing only in that the 
projection is not to SU(3) but to U(3) 

Vn,n = P r °jt/(3)[(1 ~~ a \)U n ,}i + ~r £ V n , v;fl V n+ ^ ;v V^ +fl v . ] , 
V n ,n-, v = Proj f/(3) [(l-a 2 )C/„ # + — - £ Vm^+p^;pv^+ Ap . vm ], (2-2) 

where V n ,u is the link from which the Dirac operator is to be constructed. The projection is defined 
by 

Pro]mA=A vm (23) 

which is differentiable where A is non-singular. 

The details of how the derivative is computed are described in Ref. The technique is very 
similar to the one use for stout smearing and results in similar computational cost. In principle, 
the discontinuity where the matrix A is singular could cause problems in the simulation, however, 
in practice this turns out not to be the case in the simulations presented below. Apart from the 
simulations described in these proceedings, n-HYP smearing has also been used in overlap simu- 
lations [§,0,g]. 



3. The simulation 

In order to test the performance of the smeared link, we started simulations of two degenerate 
flavors of improved Wilson fermions on 16 3 x 32 lattices. We use tadpole improved Luscher- 
Weisz gauge action. The clover coefficient is set to its tree level value csw = 1- All links of this 
operator are constructed from n-HYP links, where the coefficients are set to the standard HYP 
values, ai = 0.75, a2 = 0.6 and a 3 = 0.3. 

We generated lattices at two values of the sea quark mass. Both runs give a Sommer parameter 
ro/a = 3.8 which, using ro = 0.5fm, translates to a lattice spacing ofa« 0.13fm. One ensemble, 
labeled heavy in the following — with /3 = 7.2, K = 0.1272 — has a pseudo-scalar mass of about 
520MeV . A lighter run with j3 = 7.1, K = 0.1280 renders m PS ps 360MeV. We use multiple time 
scale Hybrid Monte Carlo[^] with Hasenbusch's mass preconditioning [jl^]. We typically ran at 
90% acceptance rate at moderate step size and unit length trajectories. 

We note in passing, that the overhead in the computation associated with the nHYP smearing — 
construction of the smeared links and the additional differentiation — is modest. For the light 16 3 x 
32 run 13% of the wall clock time is spent on this part whereas we spend two third of the time on 
the inversions and 18% on the gauge force. We believe that the cost of the nHYP smearing is more 



3 



Dynamical simulations with HYP -link Wilson fermions 



Stefan Schaefer 




Figure 1: Auto-correlation analysis of the plaquette for the light run on the 16 3 x 32, jS = 7.1, k = 0.1280 
lattice. 



than off-set by the improvement of the condition number of the fermion matrix and the reduced 
cost of the inversions. 



3.1 Auto-correlation 

In order to make sure that using smeared links does not lead to increased auto-correlation 
times, we compute T; nt for the plaquette and the plaquette constructed from smeared links which is 
less UV dominated. We use the methods of Ref . [|ll|] . For the light ensemble, we find Ti nt = 5 .2 ( 1 .4) 
and 10(4) respectively, see Fig. [T[ The heavy runs with mps ~ 520MeV have Ti nt = 6(2) and 9(3), 
which are surprisingly similar (if one ignores the large errors). However, the parameters of the 
algorithm have been tuned more carefully for the light run which might explain the absence of 
critical slowing down. 

In particular in simulations with the DBW2 gauge action, the auto-correlation time of the 
topological charge frequently are hundreds of trajectories. We test for this on a smaller lattice, 
12 3 x 24, L ps 1.4fm and mps ~ 450MeV. We measure the topological charge as defined by the 
index of the overlap operator after every fifth trajectory. The result is shown in Fig. ||. We find no 
detectable sign for increased auto-correlation. 

We take these two findings as evidence that the intuitive expectation that the smeared links do 
not have any negative effects on the auto-correlation time is indeed confirmed. 



3.2 Stability 

In order to get reliable results from a simulation one needs to be sure that it is stable. How to 
check for stability is matter of some debate. We have tried to look for signs of meta-stabilities as 
put forward by Refs. [12] and [13] and the distribution of the lowest eigenvalue of the Hermitian 
Dirac operator advocated in Refs. [14] and [15]. 
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Figure 2: Monte-Carlo history of the topological charge. Each unit on the x-axis corresponds to five trajec- 
tories of unit length by which the measurements are separated. 
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Figure 3: Monte-Carlo time history of the plaquette from a cold start. The horizontal line denotes the 
average from the hot-start. There is no sign for instabilities from a first-order phase transition 



First to the meta-stabilities: the method is to compare observables (typically the plaquette) as 
a function of Monte-Carlo time from a stream which is started from a hot configuration with a 
stream which starts from a unit, cold configuration. The result for the light run is shown in Fig. ||. 
We observe that the cold start reaches the average of the plaquette from the hot start after about 50 
trajectories and there is no sign for intermediate meta-stabilities. 

Since this is reassuring, we can now turn to the distribution of the lowest magnitude eigenvalue 
of the Hermitian Dirac operator Q = y=,D. Ref. [14] uses this distribution as the main indicator for 
the stability of a simulation. Due to the singularity in the force where Q has a zero eigenvalue, 
instabilities occur once these surfaces are crossed too often during the generation of the ensemble. 
A situation where the distribution of the lowest eigenvalue is far away from zero (where far is 
defined in units of the width of the distribution) is interpreted as stable. 

In Fig. |] we have plotted the distribution of |Ao|, the magnitude of the eigenvalue of Q with 
the smallest absolute value. It contains the result for the light and the heavy run along with its bare 
PCAC quark mass. For the heavy ensemble, we find a mean of the distribution pL = 49MeV and a 
width of a = 3MeV, the width is defined as in Ref. [14]. The light run has a median gap of 23MeV 
and also a width of 3MeV. As argued in Ref. [14] the width is independent of the quark mass. Also 
both of our runs fulfill 3 a < jJ., the criterion of stability given in that paper. If one assumes that a 
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Figure 4: Distribution of the modulus of the lowest eigenvalue of the Hermitian Dirac operator y$D for the 
two large volume runs. 
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Figure 5: Scaling of the width a of the distribution of the lowest eigenvalue of Q for various volumes and 
lattice spacings. 



is independent of the sea quark mass, and also assumes that the lower bound of stability is given 
by this criterion, one gets that a pion mass of 240MeV is reachable with our setup without running 
into problems with the stability. As a side remark, we notice that the bare quark mass is quite 
close to the median of the gap, however, some unknown renormalization constants are needed for 
a quantitative comparison. 

In the same paper [14] it has been argued that the combination aV /a should be a scaling 
quantity. Our result for a selection of our runs is displayed in Fig. |5[ We find a value of about 
0.7 whereas the original publication^] found for thin link unimproved Wilson fermions values 
around 1. However, since we use the bare width a in this plot, a meaningful comparison is again 
not possible. 



4. Summary 

At this conference, we presented new results from simulations with nHYP smeared improved 
Wilson fermions in 2fm boxes at a coarse lattice spacing of about 0.1 3fm and light pion mass 
tops ~ 360MeV. We showed that the simulations are stable in this region where ordinary thin link 
Wilson simulations are at least very difficult. Contrary to special gauge actions (e.g. DBW2) which 
suppress dislocations, smeared links have no negative impact on auto-correlation times. Given the 
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experience of quenched simulations we also expect an improved scaling behavior. This however, 
is subject of future studies. 



If one assumes the stability bound of Ref. [ |14| ] to be correct, we could simulate at pion masses 



down to 240MeV without running into problems with the stability. This would correspond to a 
mpsL « 2.4. From the scaling arguments of the same paper one would conclude that all reasonable 
parameter values in the p-regime are accessible even at such a coarse lattice spacing. 

5. Acknowledgments 

We thank the Zeuthen computer center of DESY for providing us with access to their linux 
clusters. This research was partially supported by the US Dept. of Energy. 

References 

[1] M. Albanese et al. [APE Collaboration], Phys. Lett. B 192 (1987) 163. 

[2] A. Hasenfratz, R. Hoffmann and S. Schaefer, JHEP 0705, 029 (2007) [arXiv:hep-lat/0702028]. 

[3] A. Hasenfratz, R. Hoffmann and S. Schaefer, PoS(LATTICE 2007)104 

[4] C. Morningstar and M. J. Peardon, Phys. Rev. D 69 (2004) 054501 [arXiv:hepdat/03 1 1018]. 

[5] A. Hasenfratz and F. Knechtli, Phys. Rev. D 66 (2002) 010001 [arXiv:hepdat/0 103029]. 

[6] T. DeGrand and S. Schaefer, arXiv:0708.1731 [hep-lat]. 

[7] A. Hasenfratz, R. Hoffmann and S. Schaefer, arXiv:0709.0932 [hep-lat]. 

[8] A. Hasenfratz, R. Hoffmann and S. Schaefer, PoS(LATTICE 2007)309 

[9] C. Urbach, K. Jansen, A. Shindler and U. Wenger, Comput. Phys. Commun. 174 (2006) 87 
[arXiv:hep-lat/0506011]. 

[10] M. Hasenbusch, Phys. Lett. B 519 (2001) 177 [arXiv:hep-lat/0107019]. 

[11] U. Wolff [ALPHA collaboration], Comput. Phys. Commun. 156 (2004) 143 [Erratum-ibid. 176 
(2007) 383] [arXiv:hep-lat/0306017]. 

[12] F. Farchioni et al, Eur. Phys. J. C 39, 421 (2005) [arXiv:hep-lat/0406039]. 

[13] F Farchioni et al, Eur. Phys. J. C 42 (2005) 73 [arXiv:hep-lat/0410031]. 

[14] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and N. Tantalo, JHEP 0602 (2006) 01 1 
[arXiv:hepdat/05 12021]. 

[15] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and N. Tantalo, JHEP 0702 (2007) 056 
[arXiv:hepdat/06 10059]. 



7 



